Entropy production and wave packet dynamics in the Fock space of closed chaotic 

many-body systems. 
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Highly excited many-particle states in quantum systems such as nuclei, atoms, quantum dots, 
spin systems, quantum computers etc., can be considered as "chaotic" superpositions of mean-field 
basis states (Slater determinants, products of spin or qubit states). This is due to a very high level 
density of many-body states that are easily mixed by a residual interaction between particles (quasi- 
particles). For such systems, we have derived simple analytical expressions for the time dependence 
of energy width of wave packets, as well as for the entropy, number of principal basis components 
and inverse participation ratio, and tested them in numerical experiments. It is shown that the 
energy width A(t) increases linearly and very quickly saturates. The entropy of a system increases 
quadratically, S(t) ~ t 2 at small times, and after, can grow linearly, S(t) ~ t, before the saturation. 
Correspondingly, the number of principal components determined by the entropy, N pc ~ exp{S(t)), 
or by the inverse participation ratio, increases exponentially fast before the saturation. These results 
are explained in terms of a cascade model which describes the flow of excitation in the Fock space 
of basis components. Finally, a striking phenomenon of damped oscillations in the Fock space at 
the transition to an equilibrium is discussed. 

PACS numbers: 03.67.Lx, 05.45.Mt, 24.10.Cn 
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I. INTRODUCTION 

Highly excited many-particle states in many-body sys- 
tems quite often can be presented as "chaotic" super- 
positions of shell-model basis states, see recent calcula- 



, multicharged ions H, nu- 
,|). The origin of this phe- 



tions for complex atoms 
clei || and spin systems 
nomenon relates to a very high density of many-particle 
energy levels, which increases drastically with an increase 
of energy. Indeed, the number N of combinations in 
the distribution of n particles (or quasi-particles) over m 
"orbitals" (single-particle states) is exponentially large 
(N ~ m\/nl(m — n)\ in a Fermi system). Therefore, 
the spacing D between many-body levels is exponentially 
small and a "residual" interaction V between the parti- 
cles can mix a huge number of the basis states of a mean 
field Hq (Slater determinants), when forming exact eigen- 
statcs of the total Hamiltonian H = Hq + V. 

The onset of chaos for highly excited states, as well 
as for many-particle spectra, has been recently studied 
in great detail in terms of the Two-Body Random In- 
teraction (TBRI) model which was invented about three 
decades ago In this model all two-body matrix ele- 
ments are assumed to be independent and random vari- 
ables, therefore all dynamical correlations are neglected. 
Thus, the TBRI model is essentially the random matrix 
model, however, it differs from standard Random Matrix 
Models where the two-body nature of interaction is not 
taken into account (see, e.g. |f7|-0| and review [flC 



One of the important results obtained recently [ [Tl| in 
the frame of this model is the Anderson-like transition 
which occurs in the Fock space determined by many- 
particle states of Hq. The critical value V cr for this tran- 
sition is determined by the density of states pf = dj 1 of 
those basis states which are directly coupled by a two- 
body interaction [|^,[ll|-|l3|] . When the interaction is 
very weak, Vq <C df , exact eigenstates are delta-like 
functions in the unperturbed basis, with a very small ad- 
mixture of other components which can be found by the 
standard perturbation theory. With an increase of the 
interaction, the number of principal components N pc in- 
creases and can be very large, N pc S> 1. However, if the 
interaction is still not too strong, tt~ 2 yJdfD <C Vq < df 
H , the eigenstates are sparse, with extremely large fluc- 
tuations of components. In order to have ergodic eigen- 
states which fill some energy range (see below) , one needs 
to have the perturbation large enough, Vq ^> df (for a 
large number of particles this transition is sharp and, in 
fact, one needs the weaker condition, Vq > df ). 

Above the threshold of chaos, Vq > df, the number of 
principal basis components in an eigenstate can be esti- 
mated as N pc ~ T/D where T is the spreading width of 
the strength function. In such chaotic eigenstates any ex- 
ternal weak perturbation is exponentially enhanced. The 
enhancement factor can be estimated as y/ ' N pc oc \j \J~D, 
see e.g. |l4| and references therein. This huge enhance- 
ment have been observed in numerous experiments when 
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studying parity violation effects in compound nuclei (see, 
for example, review fig]). 

In recent works ]lq , |l7| ] the theory of many-body chaos 
has been extended to quantum computers. Since in this 
case the density of energy levels is extremely high, it 
is often impossible to resolve particular many-body lev- 
els. This happens for the injection of an electron into a 
many-electron quantum dot, for the capture of an ener- 
getic particle by a nucleus or atom, or for different models 
of quantum computer with a large number of interacting 
qubits (spins). In this case the approach based on the 
study of stationary chaotic eigenstates turns out to be 
not an adequate one, and one should consider the time 
evolution of a wave function and entropy [ p~7|| . 

In contrast to the study of spectra and eigenstates, the 
analysis of the evolution of wave packets in random ma- 
trix models is based mainly on numerical results. In the 
first line, one should mention the numerical study |ll| 
of Band Random Matrices that describe quasi-lD disor- 
dered models with finite number of channels. Recently, 
the attention has been paid to the so-called Wigner Band 
Random (WBRM) model (T^] which is used in the study 
of generic properties of strength functions in dependence 
on the strength of interaction 
Refs. 



21. In particular, in 



22| , |23[ the problem of the quantum-classical cor- 
respondence for the time evolution of wave packets was 
under close study. We note that the WBRM model serves 
as convenient random matrix model for different quan- 
tum systems, and in many aspects can be compared with 
the TBRI model @. 

In this paper we study generic properties of the evolu- 
tion of wave packets in the energy shell, paying the main 
attention to the time-dependence of the entropy, width 
of packets, and inverse participation ratio. We derive 
analytical estimates for these quantities and check nu- 
merically our predictions with the use of the TBRI and 
WBRM models. 



II. TIME EVOLUTION OF CHAOTIC 
MANY-BODY STATES 



Exact ("compound") eigenstates \k) of a many-body 
Hamiltonian H — Hq + V can be expressed in terms of 
simple shell- model basis states |/) (eigenstates of Hq), or 
by products of qubits in a model of quantum computer: 



ifc) = E c / fc) i/) 



i/) = 



.a~f \" vacuum") 



(1) 



Here af are creation, or spin-raising, operators (in the 
latter case, the ground state \"vacuum") corresponds to 
the situation with all spins down), and are compo- 
nents of compound eigenstates |fc) formed by the residual 
interaction V. 



In what follows, we consider the time evolution of 
the system for the case when compound eigenstates are 
chaotic. By this term we mean that the number of prin- 
cipal components is very large, yjN pc ^> 1, and the com- 

(k) 

ponents Cj- can be treated as uncorrelated amplitudes 
with the gaussian distribution around their mean values, 
see details in Bj . Let us assume that initially (t = 0) the 
system is in a specific basis state |0) (with certain orbitals 
occupied, or, in the case of a quantum computer, when 
the state with certain spins "up" is prepared) . This state 
can be presented as a sum over exact eigenstates of the 
total Hamiltonian H: 



|o> = £ c o \k) 



(2) 



Then the time-dependent wave function reads as 

*(*) ^ E C * )C ? ] I/) ex P H£ (fe) t). (3) 
kj 

Here the sum is taken over the compound states k and 
basis states /, and we set % = 1. 

The probability W = \A Q \ 2 = | (0|*(i)) | 2 to find the 
system in the initial state is determined by the amplitude 

A Q = (0| exp{-iHt)\0) = E |Cf } I 2 exp(-iE<*>t) ~ 



J dEP {E)cxp{-iEt). (4) 



Here we replaced the summation over a large number 
of the eigenstates by the integration over their energies 
E = E^ k \ and introduced the strength function (SF) 
Pq(E) which is also known in the literature as the lo- 
cal spectral density of states, LDOS, 



P (E) = \ci k) \ME) 



(5) 



with p(E) as the density of exact eigenstates. 

As one can see, the probability Wo is entirely deter- 
mined by the strength function (||) . It is well known that 
in many applications this function has the Breit- Wigner 
(BW) form with a half- width Tq = 2irpH 2 . In our case 
p = pf is the density of directly coupled basis states and 

Hj = |-ffo/| 2 is the variance of the non-zero off-diagonal 
elements of H , defined by the residual interaction V. One 
should remind that in real situations the second moment 
A|; of the SF is always finite due to a finite range of 
interaction in the energy representation. Therefore, the 
Breit- Wigner form of the SF can occur for a finite en- 
ergy range only, determined by the energy width of the 
interaction. 

However, it was recently shown (see e.g. |8|,[l3| and ref- 
erences therein) that if r defined by the above expres- 
sion, is of the order (or larger) of the mean-square-root 
width A E of the SF itself, the form of the SF in the TBRI 
model is very close to the Gaussian. Strong deviations 
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of the SF from the BW-dependence have been observed 
numerically when studying the structure of the SF and 
eigenfunctions of the Ce atom Also, numerical data 
H have revealed that the form of the SF in nuclear shell 
models is much closer to the Gaussian rather that to the 
BW. This results from the fact that the three orbitals s, 
d 3 / 2 and d 5 / 2 included into the nuclear shell-model calcu- 
lation H , have close mean-field energies and the residual 
interaction V plays the dominant role in the Hamiltonian 
matrix. 

Recent analytical results jl3j for the TBRI model al- 
low to describe the whole transition for the SF from the 
BW regime to that of the Gaussian. This model is char- 
acterized by two-body random matrix elements which 
determine the residual interaction V between n Fermi- 
particles occupying m orbitals (single-particle states), see 
details and references, for example, in f|,Q. It was 
shown that in general case the SF is given by the fol- 
lowing approximate expression psLpjl, 



Po(E) 



1 



T(E) 



2ir (E + 6 (E) - EY + r 2 (£)/4' 



(6) 



Fermi-particles the width turns out to be the same 
for any basis state |0). 

Let us first start with the probability Wo (t) of the sys- 
tem to stay in the initial state. In two limit cases of small 
and very large times, the dependence Wo is shown p7 26 
to be of the following forms, 



and 



W (t) = exp(-A E t 2 ) 



W (t) = Ce X p(-T p t); 



t<-rf (9) 



A 2 

E 



t » TY ■ (10) 



A 2 

^E 



Here T p is the imaginary part of the pole of the SF (see 
Eq.(^J) in the complex energy plane. In the case when 
the SF has the standard BW form, we have the obvious 
relation T p = To where the latter is given by the Fermi 
golden rule. In other limit of a strong interaction, when 
the SF has the Gaussian form, the expression for T p is 
not simple. The transition from one regime of the time 
dependence of Wo(t) to another is schematically shown 
in Fig.l. 



T(E) ~ 2n\H 0f fpf(E). 



(J) 



One can see that the general expression (g) is of the Breit- 
Wigner form, however, with T(E) as some function of 
the total energy E. Here 80(E) is the correction to the 
unperturbed energy level Eq due to the residual interac- 
tion V, and Pf(E) is the density of those basis states |/) 
which are directly coupled by the interaction Hof with 
the initial state |0). It was shown jl3| that for the TBRI 
model the function T(E) has the Gaussian form with the 
variance which depends on the model parameters. In the 
case of relatively small (but non-perturbative) interac- 
tion, r = 2npfH'j <C A E , the function T(E) is very 
broad (i.e. it does not change significantly within the 
energy intervals ~ T and Ae) and can be treated as the 
constant, T(E) ~ To- In the other limit case of a strong 
interaction, To > Ae, the dependence T(E) in (^) is the 
leading one, and the slow dependence in the denomina- 
tor can be neglected. One should note that the simple 
expression (f?j) for the width T(E) has to be modified in 
this limit, see details in 



The second moment A E of the SF can be found from 



the equation A E 



E #< 2 



/¥<> 



<>/■ 



Here the summation is 



taken over the off-diagonal matrix elements Hof which 
couple the initial state |0) with the others |/). For the 
TBRI model the analytical expression for Ae has been 
obtained in Hi, 



A 



1 



— V n(n — l)(m — n)(m — n + 3) 



(8) 



where V 2 is the variance of the off-diagonal matrix ele- 
ments of the two-body residual interaction V. In fact, for 



lnW (t) 




-0.15 



FIG. 1. Schematic time dependence of Wo(t) for 
T p — 0.5, Ae = 1.2; the point t c where the dependence (g) is 
changed to (tiB is t c = T p /A% ~ 0.17. 



Now we estimate the probabilities Wf — \Af(t)\ 2 to 
find the system in other basis states. For a very small 
time we have, 



^ = |(/| e -^|0>|^|^o/| 2 t 2 



11) 



Note, that only the states directly connected to the ini- 
tial state are populated at this time scale. One can esti- 
mate the population of these states for a larger time by 
substituting the time-dependent wave function \I/(t) = 
Ao(t)ipo + J2f Af(t)ipf into the Schrodinger equation, 



dA 



m—^- = Hof (t)Ao + J2 H fk it) A k ■ 



(12) 
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Here k, f ^ and H^f(t) = Hqj exp(iu>oft). Note that 
the second term in the right-hand-side may be treated as 
a random variable with the zero mean value. Indeed, 

The variance 
I A 12 



ife oc 



H ok , therefore, H fk (t)H { 



0. 



H 



/fe(*)| 2 EfcK 1 /, 



of this term is ^ fe \Hfk(t)A}~\ 2 — 
\Hfk(t)\ 2 (l — |^4o| 2 )- Comparing this with the first term 
in the right-hand-side of Eq. ( ft2| ) , one may conclude that 
the second term is not very important for small times 
when Ao(t) ~ 1. Neglecting the second term and assum- 
ing = exp(— rt/2) which is valid for T << A E , 



we obtain 17|, 



strength function (see Eqs.(|5|, ^|)), with a slightly larger 
spreading width. 

Note that similar expression for Wq contains the term 
For gaussian fluctuations of the components 



\C\ 



(fe)|4 



C ( o>, one can get |Q KJ | 4 = 3 (\C ( >\ 2 ) 2 which is the 
known result in the Random Matrix Theory (29) . There- 
fore, if the number of principal components N pc in the SF 
is very large, the probability to find the system in the ini- 
tial state |0) at large time is, at least, three times larger 
than the probability to find the system in any other state 
|/), see Figs. 2,3. 



w f = \ H of\ 2 



\H< 



o/l 



J 0f 



r 2 /4 



|A)W|e^ 0/ *cft 



,(iwa t -Tj2)t _ l 



(13) 



where cjo/ — Ef — Eq. This approximate estimate shows 
that only the basis states within the energy interval T 
can be substantially populated (if T > Ae, this energy 
interval is equal to Ae)- 

For large time, the result is different for perturbative 
and chaotic regimes. In the perturbative regime the ex- 
pression ( [T^ ) for Wf is the final one. In the chaotic regime 
the asymptotic expression for t — * oo can be obtained in 
the following way. The projection of ^(t) (see Eq. (g)) 
onto the state / gives 



w f (t)= w s f + wf uct (t), 



(14) 



if i 



E 

k 



(*0|2|£i(fc)|2 
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FIG. 2. Asymptotic distribution Wf for the case when the 
strength function is of the Breit-Wigner form in the TBRI 
model. The parameters are: n = 6, m = 12, V 2 ~ 0.003, 
To ~ 0.50, A_b ~ 1.16, with the average over N g = 10 ma- 
trices with different realization of random two-body matrix 
elements, see in the text. 



2tt P (E - E f ) 2 + (r t /2) 2 



(15) 



0.004 



Here the result is written for the case when the SF has 
the BW form both for initial |0) and final state |/) with 
the corresponding half- widths Tg and Tq . In this case the 
resulting form of is, approximately, again the Breit- 
Wigner with the new half-width T t 
However, if Tq > Ae, the form of w s f is close to the 
Gaussian with the variance (A#) 2 



F° 
1 o 



2r n . 



The term 



fluct 



2A| 



0- 



(t) can be written in the form, 



: fluct 



(*)= E 



E (p) )t). 



At large time, t 
late and one can put Wj luct (t) 



(16) 

oo, the terms in the sum rapidly oscil- 
0. 



Thus, asymptotically 
the distribution of the components in the time-dependent 
wave function is close to that given by the form of the 



W f 



0.003 




FIG. 3. Asymptotic distribution Wf for the case when 
the strength function is close to the Gaussian. The only dif- 
ference from Fig. 2 is the interaction strength, Vq ~ 0.083, 
correspondingly, Tq ~ 10.5 and Ae ~ 5.8; the average was 
taken over N„ — 50 matrices. 
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In these figures the distribution of probabilities Wf in 
the TBRI model is shown after a very long time t — 40 for 
two different strengths of interaction (in fact, the time t 
is measured in units h/do where do — < f-s+i — e s > is the 
mean level spacing between single-particle energies e s ). 
In both cases n — 6 Fermi-particles occupy m = 12 or- 
bitals, therefore, the total number of many-particle states 
(the size of the Hamiltonian matrix) is N = 924 . The 
distribution of e s is taken random with do = 1. Two- 
body matrix elements are taken as gaussian random en- 
tries with the zero mean and variance V 2 , and in order to 
reduce the fluctuations, the distribution Wf is obtained 
with an average over a number N g of matrices with dif- 
ferent random realizations. Initially, only one basis state 
n = 462 was populated at the center of the energy spec- 
trum, in order to avoid the asymmetry of the distribution 
in the basis representation. 

The two different values of Vb, for which the distribu- 
tions Wf are obtained, are chosen in such a way that in 
one case, see Fig. 2, the strength function has the Brcit- 
Wigner form, and in the other, the form is very close to 
the Gaussian (Fig. 3). We should remind that the above 
two forms occur in the energy representation, however, 
the results are shown in the basis representation. These 
two representations are related through the density of 
states which is known to be of the Gaussian form for 
large number of particles and orbitals Bp^l ■ 

One can see that in both cases, Fig. 2-3, the probability 
to stay in the original basis state is much larger than in 
the nearest ones. Comparing with the result of the stan- 
dard random matrix theory, one can say that there is a 
noticeable difference (namely, the enhancement factor in 
Fig. 3 is about 2.3, instead of 3.0). 



III. THE CASCADE MODEL 

One of the important question is how the entropy of 
quantum isolated systems increases in time at the tran- 
sition to equilibrium. It is natural to define the entropy 
of a many-body state through the Shannon entropy, 



S(t) 



22 w k In Wk 



-WolnVFo - 

Wf In wf 



(17) 



Here Wo(t) = |^4o(i)| is the probability for the system 
to be in the initial state, and Wf (t) = \Af (t)\ is the 
probability to be in the basis state |/). In what follows 
we assume that the initial conditions are Wo(0) = 1 and 
Wf(0) = 0, therefore, the entropy is equal to zero for 
t = 0. 

In order to study the evolution of a many-particle sys- 
tem with two-body interaction, it is convenient to intro- 
duce sub-classes for all basis states in the following way. 
The first class contains those Ni basis states which are 



directly coupled with an initial state by the two-body 
interaction given by matrix elements Hof of the inter- 
action. Correspondingly, the second class consists of N2 
basis states which are coupled with the initial one in the 
second order of the perturbation, this coupling is deter- 
mined by Hf) a H a f , etc. 

Let us first consider the evolution for a large time 
t ^> (below we assume the BW shape of the SF). 

For this case the probabilities of the states in different 
classes can be determined by the "probability conserva- 
tion equations" , 



dW 
dt 



-YW a 



dWi 
dt 



rw - 1W1 



dWk 
dt 



= TW k -i - TW k 



(18) 



Here Wk is the probability for the systems to be in the 
class k. The first term TWk-i in the right-hand-side of 
( |l8| ) is responsible for the flux from the previous class, 
and the second term TWk describes the decay of the 
states in the class k, into the next class fc+1. We assume 
that the probability of the return to the previous class 
can be neglected. This is a valid approximation if the 
number of states N^+i in the next class is large in com- 
parison with Nu of the previous class. This approach can 
be compared with those based on the Caley tree model 
]ll| where the flow from each state goes into other M 
states, therefore, N k ~ M k with M > 1. Note that here 
we consider a system which is far from the equilibrium. 
Indeed, if the system is in the equilibrium, the probabil- 
ities of all states within the energy shell defined by the 
relation \Ef — Eq\ < T, are of the same order, Wf ~ N~^, 
with N pc as the total number of states inside the energy 
shell. Therefore, in order to neglect the return flux, one 
needs the condition Wf = Wk/Nk ^> 1/N pc to be fulfilled. 
Equations (Il8[) have the simple solution, 



W = exp(-rt) 



Wn = { -^exp(-Tt) = ^Wo. 



(19) 



The maximal probability W„ — ^yexp(— n) as l/v2~7ra 
to be in the class n determined by the condition ^f - — 
, occurs for t = n/T , therefore, this solution ( J19| ) can 
be considered as a cascade in the population of different 
classes . Indeed, at small times t <C t = 1/T the system 
is practically in the initial state, at times t w r the flow 
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spreads into the first class, for t = nr it spreads into the 
n— th class, etc. For an infinite chain one can easily check 
the normalization condition. 



n=0 



W n = exp(-rt) 



oc 

E 

n=0 



(Tt) n 



= 1. 



(20) 



IV. TIME DEPENDENCE OF THE ENTROPY 

The above expressions allow us to find the time- 
dependence of the entropy, 



S(t) 



n=0 



Wn 



TthxM +Tt- e 



-rt 



y, (rt) n 1 (rt) n 



(21) 



n=0 



where w/ ~ stands for the population of basis states 
of the class n with N n as the number of states in this 
class (in fact , for t ~ rvr one needs to count only the 
states inside the energy shell since the population of the 
states outside the energy interval with \Ef — E \ > T, 
is small, see Eq.(|l3|)). Here we have used the relations 

N n = M n and ^7T" n = Ttexp(Tt). Two last terms 

n=0 

in the right-hand-side of Eq. (|l| ) turn out to be smaller 
than the first one, therefore, one can write, 



S(t) S3 K In M[l + /(*)] 



(22) 



with some function f(t) <C 1 which slowly depends on 
time. 

In this estimate for the increase of entropy, we did 
not take into account the influence of fluctuations of 
Wf. One can show that for gaussian fluctuations of the 
coefficients Af with the variance given by their mean- 
square values, for large number of principal components 
N pc (t) = exp(S(t)) the entropy should be corrected by a 
small factor of the order of In 2 (see, for example, |3C|). 

If one neglects the second term in ( 22 ) , we obtain a lin- 
ear increase of the entropy, which means that the number 
of principal components N pc (t) increases exponentially 
fast with time. This behavior can be compared with 
a linear increase of dynamical entropy S c i (t) in classical 
chaotic systems where S c i (t) was found to be related to an 
exponential divergence of close trajectories in the phase 
space (S c i(t) oc Xt with A as the Lyapunov exponent, see, 
for example pij ). The non-trivial point is that the linear 
increase of entropy also occurs for systems without the 
classical limit, see recent paper B3]. 

Note, that for a small time the function Wo (t) has the 
form Wo(t) = exp(— A^t 2 ) (see ([|)), not the exponen- 
tial dependence exp(— Tt) . Therefore, one should modify 
the expression for the entropy in order to make it valid 



for small times. For this, we replace Tt in Eq. (]Ts|) by 
a more accurate expression, — ln(W'o), which gives A|. t 2 
for small time t <§; T/A% and Tt for large time t ^> T/A 2 E 
, therefore, 



ni 



(23) 



It is easy to check that the normalization condition is 
fulfilled again, 

E W n = E , W = WoexpiHW- 1 )) = 1. 



n— 1 n— 1 



For the entropy one obtains 



n=l 



Wn 



(24) 



(25) 



At small time Aj^ the entropy is given by two 

terms, n = and n = 1 (direct transitions), therefore, 



S{t) = -W Q {t)\nWa{t) - Wi{t)\n 



Alt' 



Wi 

iv7 



A 2 E t 2 1 + In ' 



(26) 



This expression can be compared with the direct calcu- 
lation based on the relation Wf — |iJo/^| 2 , 



S(t) = A 2 E t 2 +t 2 J2\H af \ 2 \n 



1 



H 2 t 2 



(27) 



There is an agreement between Eqs.(|26|) and (|27j) since 

Ai = E/i^o/i 2 =^ii^;. 

Let us now discuss the whole time dependence of the 
entropy, including large times t when the system is close 
to the equilibrium. In a finite system of particles any ba- 
sis state can be reached, starting from an initial state, 
in several "interaction steps" {HQ a H a pH/3-y . . .) . For 
example, in the system of n = 6 particles three steps 
is needed since the two-body interaction can not move 
more than two particles from one basis state to another. 
If the number of classes n c is finite, the states in the 
last class do not decay (there is no term —TW nc in the 
last equation in flis|)), and the probability to be in the 
last class is determined from the normalization condition 
W Ua = 1 — Y^k'=Q 1 Wk- The additional condition is that 
the considered basis states should be inside the energy in- 
terval, \Ef — Eq\ < min(T, AE) , thus limiting the num- 
ber N n of the basis states in each class. Note that the 
value of W n is restricted from below by the equilibrium 
relation W n > N n /N pc . These limitations make "exact" 
expression for the entropy very complicated. Instead, 
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we can propose the following simple expression which is 
approximately valid in systems with a small number of 
classes (n c ~ 1), 

S(t) = -W (t)lnW (t) - (1 - Wb(*))ln ( ^77^ 

(28) 

This expression takes into account the normalization con- 
dition ^2 Wf — 1 — Wo and has a reasonable behavior 
J¥o 

for both small and large times. 



V. NUMERICAL RESULTS FOR THE ENTROPY 

Now we compare the obtained analytical expressions 
with numerical data for the TBRI model. For the case 
when the strength function (SF) has the Breit-Wigncr 
form, the time dependence of the entropy is shown in 
Fig. 4 for the parameters of Fig. 2, n = 6, m = 12, Vq 2 s» 
0.003, To ~ 0.50, As sw 1.16, with the average over 
N g = 2 Hamiltonian matrices. 



S(t) 




FIG. 4. Entropy versus time for the TBRI model in the 
case when the SF has the standard Breit-Wigner form. The 
parameters are the same as in Fig. 2, with N g = 2. The cir- 
cles stand for numerical data, the solid curve is the analytical 
expression (p8|), and the dashed line represents a linear slope 
according to an approximate expression ( fjo|) . In the inset the 
same is shown for a smaller time scale. 



The number M of the basis states directly coupled by 
the random two-body interaction, is determined by the 
expression 0, 

n(n — 1)(to — n)(m — n — 1) 



M = n(m — n) + 



(29) 



where the first term gives the number of one-particle 
transitions, and the second stands for two-particles tran- 
sitions. In our case of n = 6 particles and to = 12 or- 
bitals, the total number of basis states is N = 924 and 



M = 261. The effective number of classes in the cascade 
model can be determined from the relation M" c = N. 
This gives n c = In N/ In M w 1.2. Thus, we can use 
the simple expression ( p8| ) to describe the dependence of 
the entropy on time analytically. The data in Figs. 4-5 
demonstrate an excellent agreement between the numer- 
ical and analytical results. 

We should note that the theoretical dependence ( psj ) 
which gives quite good approximate description of the 
data on the whole time scale, has the parameter N pc (ef- 
fective number of principal components in the stationary 
distribution Wf(t — > 00)) which is related to the limiting 
value of the entropy, N pc = In 5(00). It can be estimated 
analytically as discussed above from the width of the en- 
ergy shell; when plotting the solid curves in Figs. 4-5 we 
have used the exact value found numerically. 

To avoid confusion, we should explain that the actual 
number of classes in the case of n — 6 particles and 
m = 12 orbitals is equal to 3 since all basis states can be 
populated in the third order in the two-body interaction. 
However, the number of states in the second, k = 2, and 
third, k = 3, classes are much smaller than it may follow 
from the exponential relation Nk = M k (in practice, this 
relation may be justified for a large number of particles 
only). This is the reason why the one-class formula fl2S| ) 
works so well. 

It is also instructive to compare the entropy with the 
linear time dependence, 



S(t) = Tt InM 



(30) 



that stems from Eq.(p2|) if the first term is taken only. 
This dependence globally corresponds to the data on 
some time scale, however, the actual dependence for S(t) 
clearly differs from the linear one (see inset in Fig. 4). 
Note that the quadratic increase of energy occurs on a 
very small time scale only. As for the oscillations of the 
entropy for a very large time close to an equilibrium, this 
phenomenon will be discussed below. 




FIG. 5. Time dependence of the entropy for the TBRI 
model when the strength function is of the Gaussian form: 
(a) n = 6, m = 12, Vo « 0.083, To w 10.5, A E » 5.8, and (b) 
n = 7,m = 13, V w 0.12, To w 14.6, A E « 8.13. Circles are 
numerical data for N a — 2, solid curves stand for the approx- 
imate expression j28|), and dashed lines represent the linear 
dependence (|3l|). 
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For a strong interaction, when the form of the SF is 
very close to the Gaussian, numerical data are reported 
in Fig. 5 for n = 6, m = 12 and for n = 7, m — 13. The 
interaction strength is chosen in order to have the same 
ratio r /A_B sa 1.8, as in Fig. 4. 

In this case the half-width of the strength function is 
determined by since Tq — 2-n:pfH^ is larger than 



Ae- As a result, the role of T in the expressions (|22j) and 
( p8| ) plays the width Ae- In both cases numerical data 
give strong evidence of a linear increase of the entropy, 



S(t) = Ae t In M, 



(31) 



before the saturation. It is clearly seen that this analyt- 
ical estimate gives a correct value for the slope of S(t). 
The shift is due to the initial time scale where the time 
dependence is quadratic, this fact is neglected in the es- 
timate. 

It should be pointed out that the linear dependence of 
S(t) in Fig. 5 is much more pronounced than in the BW- 
rcgion (compare with Fig. 4 for which the SF is of the 
Breit-Wigner form). Our results indicate a clear differ- 
ence between the two cases related to the Breit-Wigner 
and Gaussian forms of the SF. This point is somehow 
supported by recent studies |^2|,^3| where it was shown 
that for a relatively weak interaction (ro is small in com- 
parison with Ae) resulting in the Breit-Wigner form of 
the SF, there is no detailed quantum-classical correspon- 
dence for the evolution of wave packets in the energy 
space. On the other hand, in the Gaussian region (with 
the Gaussian form for the SF), the detailed quantum- 
classical correspondence is possible [ p2"|]23| ]. In the latter 
case one can expect a linear growth for the entropy, as 
it was found in classical models J3l[]. The principal dif- 
ference between these two cases (in what concerns the 
quantum-classical correspondence and the possibility of 
the localization in the energy shell) for the first time has 
been discussed in Ref. KG]. 



the matrix elements Hq / , to the final width of the wave 
packet in the basis representation, which is defined by 
the width of the energy shell (approximately equal to 
\/2AE) . Therefore, the leading time dependence is given 
by the term 1 - W (t). 

For relatively small times before the saturation, we can 
use the simple estimate Wf ~ \Hoft\ 2 , which results in 
the following quadratic dependence, 



A 2 (i) 



/ 



no) 



n 2 

H 0f 



t 2 V Q 2 Al 



(33) 



Here Ao is some constant related to an effective band- 
width of the Hamiltonian matrix. The linear dependence 
for the width of packet, A(t) = t Vq A , corresponds to a 
generic ballistic-like behavior of wave packets found for 
the WBRM model |§. 

Note that the band-width of the Hamiltonian matrix 
can be much larger than the final width of the wave 
packet due to the dependence of the latter on the in- 
teraction strength. Therefore, the linear increase of A(t) 
can be very fast, and quickly it becomes saturated on the 
time scale of the applicability of the expansion in Ht, see 
Eq. @. 

Before comparing the obtained expressions with nu- 
merical data, let us first analyze the time dependence 
of the number of principal components N pc (t) for wave 
packets in the basis representation. It is natural to de- 
fine N pc through the entropy, N pc {t) — exp(S(t)). This 
definition has been widely used in different applications, 
see for example, p0| . 

The number of principal components can be also de- 
fined through the inverse participation ratio li pr : 



k 



' o 2 l2 ' 

k 



l\2k 



( In IV) 
k\k\N k 



(34) 



VI. WIDTH OF PACKETS AND INVERSE 
PARTICIPATION RATIO 

The width of the wave function in the basis represen- 
tation can be measured through the variance, 

A 2 (t) = £ (n f n f \A f (ttf = (1 - W (t)) \A f (t)\ 2 ; 
/ 



ia/ (t)r 



(n -n f f\A f (t)\' 

Et \A f (t)\ 2 



(32) 



where nj and uq label corresponding basis states, 
/ 7^ 0, and we have used the normalization condition 
E/ \ A f Ml 2 = 1 - W (t) . The function |A/ (t)| 2 is a 
slow function of time, it changes from the effective band- 
width of the Hamiltonian matrix, entirely determined by 



Here we used Eq.(p3|) for Wk- The sum in Eq.(g^) gives 
the following result for the infinite number of classes (this 
may be a reasonable approximation for time t « n c r) 

(tWl^^i^), (35) 

where Iq(z) stands for the modified Bessel function and 
we used the relation N n = M n . 

This expression has the following asymptotics, 



(Mt)r 1 = in 2 (i + ^(in^v 1 ) 2 ) 



1 - 2< 2 (A.E) 2 (1 



(36) 



for small time, and 
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(Jipr (^) ) 



exp 



-2rn 



'M 



(37) 



for large time. Therefore, N pc defined through the in- 
verse participation ratio U pr also may have an interval 
of exponential increase in time (if the number of classes 
n c is not small and we can extend the summation over 
k in Eq. (|34| ) to infinity) . Here we again neglected the 
fluctuations of A f (t) which may increase the value of 
up to the factor 3. 

For a system with small number of classes on can sug- 
gest the following approximate expression, 



(V(*)) =wg + 



(1 - W f 

lipr (p®J 



(38) 



This expression takes into account the normalization con- 
dition Y] Wf = 1 — Wo and has a reasonable behavior 

for both small and large time. 

In any system with finite number of particles the en- 
ergy shell contains finite number of basis states. For a 
stationary chaotic state the number of principal compo- 
nents is estimated as N p l ~ T/D where D is the mean 
energy distance between many-particle levels (we assume 
here that the spreading width T of exact eigenstates is 
less than Ae)- In the non-stationary problem this leads 
to the saturation of N pc (t) to the value N pc (oo) ~ 27Vp*, 
and to the maximal value 5* « In N pc (oo) for the entropy 
(see above and in Ref. p7|). 

Numerical data for the TBRI model for the case when 
the SF is of the Breit-Wigner form are summarized in 
Fig. 6 for the same parameters as in Fig. 2 and Fig. 4. 
Three quantities are plotted here: the width A(t), the 
number of principal components N pc (t) = exp(S(t)) and 
Uprit) = (J2 \A f (t) I 4 )" 1 determined by the inverse par- 
ticipation ratio. For N pc (t) two curves are given, one 
is due to the analytical expression (^8|) for the entropy, 
and another is computed directly from the evolution of 
the TBRI model, with additional average over N g = 2 
number of realization of random Hamiltonian. 

From the reported data one can see that the time de- 
pendence of the width A(t) of packets is quite simple. 
Namely, on the first very short time scale the increase of 
the width is linear in time (see also inset in Fig. 6 where 
A(i) is shown on this time scale), and after, the width 
quickly saturates. This behavior is in the correspondence 
with analytical estimates discussed above, and with nu- 
merical results found in the WBRM model p2[ . 

In contrast to this time dependence, the increase of the 
number of principal components N pc (t) is very different, 
both for N pc (t) defined by the entropy and for the h pr (t) 
determined by the inverse participation ratio. Indeed, 
both N pc (t) and li pr (t) increase slowly in time before the 
saturation to their limit values. The absolute difference 
between these two quantities is not important since the 
definition of U pr is given up to some factor which is sen- 
sitive to the type of fluctuations in Af (t) . As we already 



noted, the gaussian fluctuations decrease the value of li pr 
by the factor 3. 




FIG. 6. Time dependence of different quantities for the 
TBRI model in the case of the BW-form of strength functions. 
The parameters are the same as in Fig. 4, with N g = 2. Circles 
stand for numerical data for N pc = exp(S(t)) with S(t) taken 
from Fig. 4, dots correspond to the analytical expression for 
N pc (t) with S(t) from Eq. (jiij) , triangles represent numerical 
data for the width A(t), and squares are numerical results for 
li pr (t)- The width A(t) on a smaller time scale is shown in 
the inset. 



One of the most interesting facts which can be drawn 
from these data is a big difference for characteristic time 
scales which correspond to the saturation. Indeed, if the 
width A(i) completely saturates at time t ~ 0.5, both 
N pc and li pr manifest a very slow saturation by the time 
t ~ 20. This means that the mechanism for the width in- 
crease is different from that responsible for the increase of 
the number of principal components in the wave packet. 

To explain this phenomenon let us consider the ini- 
tial time scale for the time dependence of N pc and U pri 
where one can detect an approximate linear increase of 
the entropy S(t). The very point is that at small times 
t < t the wave function has a large number of holes since 
only directly connected basis states are populated. For 
a system with a large number n of particles the fraction 
of these states is exponentially small (~ exp(— n)) due a 
two-body nature of the interaction. With an increase of 



time, for t > r 



the states in other classes start 



to be filled and the holes begin disappear. This stage 
for t < n c r corresponds to an exponential increase of the 
number of principal components and linear increase of 
the entropy. 

It is instructive to analyze the evolution of wave pack- 
ets on a smaller time scale of the ballistic spread, see 
Fig. 7. These data confirm theoretical expectations ac- 
cording to which at small times only those basis states 
that belong to the first class, are involved in the dynam- 
ics. Indeed, large gaps are clearly seen in the distribution 
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of Wf, which persist during all the time of the ballistic- 
like spread. These gaps reduce the number of the princi- 
pal components, however, they are not important for the 
calculation of A(t), see Eq.(|3^). 

An important peculiarity of the wave dynamics is that 
initially all basis states of the first class are excited. One 
can see that for a very small time t — 0.2 the whole avail- 
able region 1 < / < N is filled with approximately the 
same amplitudes Wf = \H ft\ 2 . With an increase of time, 
the amplitudes grow and form the envelope of the packet, 
in accordance with Eq.(|l5|). One should stress that the 
quadratic time dependence for the second moment of a 
packet on this short-time scale occurs not due to a linear 
spread of the front of the wave packet, but due to a spe- 
cific growth of the amplitudes of those basis states which 
are located inside the energy shell. 
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FIG. 7. Wave packet w f (t) for the TBRI model at dif- 
ferent times t — 0.2, 0.4, 0.6.0.8 for the parameters of Fig. 7. 
One particular Hamiltonian matrix is used without any addi- 
tional average. 



The remarkable effect is a kind of oscillations for all 
quantities of Fig. 6. Similar oscillations (however, not so 
strong) are also present in the time dependence of the en- 
tropy S(t), see Fig. 4. Since the period T ss 6.5 of these 
oscillations is much larger than the time scale t w 0.5 of 
the ballistic spread of wave packets, it is clear that this 
effect is entirely related to the dynamics in the Fock space 
formed by different classes. The origin of these oscilla- 
tions can be explained in terms of the cascade model dis- 
cussed in Section III. Indeed, one can expect a strong ef- 
fect of a reflection due to the finiteness of the Fock space. 
The first reflection occurs for to « n c (Fo) _1 , therefore, 
the period of oscillation is T w 2£o- One can see that 
this estimate gives the correct result for T with n c ~ 1.5. 
The latter value is close to our rough estimate n c ~ 1.2 
for an effective number of classes. 

To compare the data for l ipr defined by the inverse 



participation ratio with the analytical expression (38), 
we have plotted separately both results for a larger time 
scale t < 40 in Fig. 8. One can see that our estimate (pSSj) 



gives quite accurate description of the data on a large 
time scale up to t « 20. After this time, the saturation 
occurs and all local time dependence may be treated as 
fluctuations around the limiting value. The data pre- 
sented in this figure give strong evidence of the effective- 
ness of our analytical approach. 




FIG. 8. Time dependence of h pr (t) for the parameters 
of Fig. 6 on a larger time scale for one Hamiltonian matrix. 
Solid curve stand for numerical data, and the dashed curve 
corresponds to the analytical expression (|§8|). 

Let us now come to the case when the form of the 
strength function is close to the Gaussian. In Fig. 9 one 
can see that the analytical expressions connecting the 
probability Wo(t) to stay in an initial basis state, with 
the time dependence of the number of principal compo- 
nents N pc (t) = exp(5(i)) and k pr (t), give correct global 
description of numerical data. When comparing with the 
previous case of the BW-form of the SF, see Figs. 6, 8, 
one can conclude that for a strong interaction (when the 
Gaussian dependence for the SF emerges) the behavior 
of all quantities do not reveal strong oscillations. This 
effect is related to the fact that the time scale for the 
saturation of the width of packets is of the same order as 
for the saturation of both N pc and li pr . In such a case 
the effect of reflection in the Fock space is suppressed by 
a strong spread of packets in the energy shell. 




FIG. 9. The same quantities as in Fig. 6 for a strong 
interaction with the Gaussian form of the SF. The pa- 
rameters for the TBRI model are the same as in Fig. 5a: 
n = 6,m = 12, V ~ 0.083, r « 10.5, As » 5.8. Solid 
curves correspond to numerical data, dashed curve stand for 
the corresponding analytical expressions: Eq.(^) for U vr {t), 
and Eq.(|2S|) for S(t) in the definition N pc (t) = exp(S(t)). 
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One can see, that there is a quite strong difference in 
the time-dependence of all discussed above quantities for 
the two extreme cases of the Breit- Wigner and Gaussian 
forms of the strength function. In the BW-case the two 
effects (ballistic-like spread of packets and the cascade- 
like evolution in the Fock space) have very different time 
scales, and both these effects can be distinguished in the 
dynamics. Contrary, in the second case of the Gaussian 
form of the SF, two time scales are comparable. There- 
fore, the two effects coexist on the same time scale and, 
as a result, the global time-dependence turns out to be 
much simpler. 

VII. COMPARISON WITH WIGNER BAND 
RANDOM MATRICES 

In this Section we discuss numerical results obtained 
for the WBRM model for the same quantities as the con- 
sidered above. The dynamics of wave packets in this 
model has been recently studied in (2^] in the connec- 
tion with the problem of quantum-classical correspon- 
dence. Here, instead, we concentrate our attention on 
the correspondence between the evolution of packets in 
the WBRM and TBRI models. 

As was mentioned, the WBRM model is quite close to 
the TBRI model. It consists of two parts, one of which 
is a diagonal matrix with increasing "energies" ej and 
another is a band matrix Vij, 

H i j = e j 8 ij +V ij , (39) 

where Sij is the delta-function. In original papers JToj ] 
the "unperturbed spectrum" was taken in the form of 
the "picked fence", Cj = jD, where D = p^ 1 is the spac- 
ing between two close energies and j is a running integer 
number. We consider here the case with random values 
6j with the mean spacing D, reordered in an increasing 
way, e^+i > tj. As for the off-diagonal matrix elements 
Vij, they are assumed to be gaussian random and inde- 
pendent variables inside the band \i — j\ < b, with the 
zero mean, < Vij = >, and the variance < V? >= V 2 . 
Outside the band, the matrix elements vanish. Thus, the 
control parameters of this model are the ratio Vq/D of 
a typical matrix element to the mean level spacing, and 
the band-width b. As one can see, the first term in ( |39| ) 
corresponds to a "mean field" H$, and the interaction V 
has a finite energy range. 

For the first time the SF for the WBRM has been an- 
alyzed in Ref. It was analytically found that the 
form of the strength function essentially depends on one 
parameter q — Po ^° only. Wigner has proved that 
for a relatively strong perturbation, Vo 3> D, in the limit 
q <C 1 the form of the LDOS is the Lorentzian, 

Wbw{E) = ± T ™ , E = E — D j, (40) 

^ n -Cj T Tl rw 



which nowadays is known as the Breit- Wigner (BW) de- 
pendence. Here the energy E refers to the center of the 
distribution. The half- width Tbw of the distribution ( |40| ) 
is given by the Fermi golden rule, 

T BW = 2n Po V 2 (41) 

In other limit q 1 the influence of the unperturbed 
part Ho can be neglected and the shape of the SF tends 
to the shape of density of the states of the band random 
matrix V, which is known to be the semicircle. 

Recently, Wigner's results have been extended in pl| ] 
to matrices H with a general form of V, when the 
variance of the off-diagonal matrix elements decreases 
smoothly with the distance r = — from the prin- 
cipal diagonal. In this case the effective band size b is 
defined by the second moment of the envelope function 
f(r) . Another important generalization of the WBRM 
studied in |PH , is an additional sparsity of the matrix V, 
which may mimic realistic Hamiltonians. In such a form, 
the WBRM model is closer to the TBRI model, however, 
in the latter the sparsity of the interaction is due to a 
two-body nature of the interaction. As a result, the po- 
sitions of zero elements are not completely random as in 
the WBRM model, see details in [@,§. 

Random matrix models of the type (|3^) are very useful 
for the understanding generic properties of the SF. The 
condition for the SF to be of the BW form in the WBRM 
model has simple form plf , 

D < T BW < A fc ; A b = bD. (42) 

The left part of this relation indicates the non- 
perturbative character of an interaction, according to 
which many of unperturbed basis states are strongly cou- 
pled by an interaction. On the other hand, the interac- 
tion should not be very strong, namely, the width Tbw 
determined by Eq.(f4l]), has to be less than the width At, 
of the interaction in the energy representation. The lat- 
ter condition is generic for systems with finite range of 
the interaction V. One should stress that, strictly speak- 
ing, the BW form (|^) is not correct since its second mo- 
ment diverges. As was shown in Jl^,|l|, outside the energy 
range \E\ > Afc the SF in the model ( ]39| ) decreases with 
the energy faster than a pure exponent. 

Note, that in the TBRI model the energy scale At, is 
irrelevant since there is no sharp border of the interaction 
and At, is of the order of the whole energy spectrum. For 
this reason, instead of A^ it is more convenient to use 
the variance A^ of the SF which may have the classi- 
cal limit pb] , |22| . The latter quantity can be expressed 
through the off-diagonal matrix elements of the interac- 
tion, A| = J^j v ?j for i¥"h therefore, A| = 2bV 2 . As 
a result, we have A& = ttA e /Tbw and Eq.(^2| ) can be 
written as 
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D < T BW < Ab-v/tt- (43) 

Numerical data Q for the WBRM model show that 
on the border Ybw ~ 2A^ the form of the SF is quite 
close to the Gaussian, and this transition from the BW 
dependence to the Gaussian-like turns out to be quite 
sharp. Although the extreme limit of a very strong in- 
teraction, g > 1, (or, the same, Tbw Ae) has been 
studied by Wigner in the WBRM model (|3^) , the semi- 
circle form of the SF seems to be unphysical. Indeed, this 
form is originated from the semicircle dependence of the 
total density py{E) defined by V only, therefore, when 
neglecting the term H . Thus, the case V ^> H in terms 
of the TBRI model means that the residual interaction 
is much stronger than the mean-field part Ho, which is 
physically irrelevant. 
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FIG. 10. Time dependence of the entropy S(t) (in 
the inset), the width of the packet A(t) and the number 
of principal components h pr (t) for the WBRM model for 
the case of the BW-form of the SF. The parameters are: 
N = 924, b = 110, D = 1.0, Vo = 1.0, correspondingly, 
Tbw ~ 6.28, A E « 14.8. Solid curves for S(t) and l ipr (t) 
are numerical data, dashed curves represent analytical expres- 
sions ( p£j| ) (for S(t)) and (0) (for li pr ). Diamonds, connected 
by the solid line, correspond to numerical data for A(t). 

Numerical data for the WBRM model in the case of 
the BW-dependence of the SF are given in Fig. 10. When 
comparing with the corresponding quantities discussed 
above for the TBRI model (see Fig. 6) we should note the 
following. First of all, one can see that the simple ana- 
lytical expression (^ ) gives a correct description of the 
increase and saturation of the entropy S(t). The same 
we can say about the expression ( |38| ) for the number of 
principal components defined through the inverse partic- 
ipation ratio. 

Second, we would like to stress that the global time de- 
pendence for all quantities is quite similar to that found 
for the one-class variant of the TBRI model. Relatively 



simple structure of the Wigner band random matrices 
allows one to perform a detailed comparison of the data 
with analytical estimates. Indeed, the application of the 
relation (|3^) for the WBRM model gives, 

A 2 (t) = -t 2 V 2 b\ (44) 
3 

see also Ref. jig]. Therefore, for the parameters of Fig. 10 
we have A(i) = Bt with B = J%b 3 / 2 V » 950 which is 
in a good agreement with numerical data. We can also 
find the critical time after which the ballistic spread 
of the packet terminates. For this, we estimate the max- 
imal width A m of packet via the width A^ of the SF, 
A m w >/2Ae- This leads to the estimate t m w v6/6, 
therefore, for Fig. 3 we have t m w 21 which perfectly cor- 
responds to the data. The latter estimates have been also 
checked for other values of Vo and b, with the same good 
correspondence between simple estimates and numerical 
data. 

Comparing global time dependence of the quantities 
presented in Fig. 10 with the results for the TBRI model, 
see Fig. 6, one can see that the main difference is the type 
of oscillations for the width of packets A(t) . Namely, in 
contrast to the TBRI model where the period of oscilla- 
tions is much larger than the time scale t m of the ballistic 
spread, in Fig. 10 the period T is just defined by the bal- 
listic spread, T » 2t m . This very fact demonstrate the 
principal difference between the two models. 

Indeed, for the WBRM model there is no specific evolu- 
tion in the Fock space which is due to a two-body nature 
of interaction. Formally, the cascade model can be ap- 
plied to the WBRM model with the number of classes 
n c = 1, since all states within the energy band A& = bD 
start to be involved in the dynamics immediately. This 
means that in contrast to the TBRI model, in the WBRM 
model there is only one mechanism for the oscillations, 
namely, the reflection inside the energy shell which is 
populated ergodically, when time is running. No oscil- 
lations are detected for the number of principal compo- 
nents (the data for larger times are not shown), this con- 
firms our conclusion about one kind of the reflection from 
the edges of the energy shell. It is interesting to note 
that the number of principal components does not reveal 
noticeable oscillations on the time scale of the ballistic 
spread since on this scale the value of N pc is very small. 

The difference between these two models can be also 
seen when comparing the structure of wave packets in 
the basis representation at some instant times before the 
saturation (compare Fig. 11 with Fig. 7). In contrast to 
Fig. 7 where many "holes" can be realized in the distri- 
bution Wf, for the WBRM model the filling of the avail- 
able energy range of size 2bD occurs ergodically. In both 
cases very strong fluctuations are present, which are ex- 
pected to be gaussian, see discussion in It is im- 
portant to stress, that in order to reveal this difference, 
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we should avoid the ensemble average which washes out 
the presence of holes (if different Hamiltonian matrices 
have different unperturbed spectrum). This fact reflects 
one of basic peculiarities of the TBRI model, namely, the 
non-ergodic character of the matrices (the average over 
the spectrum inside one (very big) matrix may give com- 
pletely different result from that obtained by an ensemble 
average, see references in the review ]29|). 
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from the inverse participation ratio, numerical data man- 
ifest the oscillations, as well, with the same period as for 
the width of packets. In average, the numerical data for 
S(t) and li pr are well described by the simple analytical 
expressions relating these quantities to the probability 
Wo(t) to stay in the initial state. 

The most interesting result which can be drawn from 
the numerical data presented in Fig. 10 and Fig. 12, is a 
clear difference in the time dependence of the entropy 
S(t). Comparing the data in the insets, one can con- 
clude that linear increase of the entropy occurs for the 
case when the form of the SF is close to the Gaussian, 
but not for the case of the BW-form. Indeed, it is hard 
to indicate clear time scale of a linear increase of S(t) in 
Fig. 10, in contrast to Fig. 12 where the linear time depen- 
dence is clearly seen (apart from a very small time scale). 
This very fact may be of quite generic, since in the TBRI 
model we also see a non-linear character of the entropy 
increase in the BW-regime. 
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Fig. 10. One particular band random matrix is used without 
any additional average. 



This work was supported by the Australian Research 
Council. One of us (F.M.I.) gratefully acknowledges the 
support by CONACyT (Mexico) Grant No. 34668-E. 




0.0 0.1 0.2 0.3 0.4 0.5 

t 



FIG. 12. The same as in Fig. 10 for the case when the 
form of the SF is close to the Gaussian. The parameters 
are: N = 924, b = 110, D = 1.0, Vo = 3.0, correspondingly, 
Ybw » 56.5, A E ^45.5. 

Finally, we present the data for the WBRM in the case 
of the Gaussian form of the strength function, see Fig. 12. 
Here we also can see oscillations for the width of pack- 
ets. As for the number of principal components found 
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